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Abstract. Small-angle correlations of pairs of protons emitted in central collisions of Ca -I- Ca, Ru -|- 
■ Ru and An -I- Au at beam energies from 400 to 1500 MeV per nucleon are investigated with the FOPI 

$H ' detector system at SIS/GSI Darmstadt. Dependences on system size and beam energy are presented which 

extend the experimental data basis of pp correlations in the SIS energy range substantially. The size of the 
proton-emitting source is estimated by comparing the experimental data with the output of a final-state 
interaction model which utilizes either static Gaussian sources or the one-body phase-space distribution of 
protons provided by the BUU transport approach. The trends in the experimental data, i.e. system-size 
and beam energy dependences, are well reproduced by this hybrid model. However, the pp correlation 
function is found rather insensitive to the stiffness of the equation of state entering the transport model 
calculations. 



PACS. 25.70.Pq, 25.75.Gz 



1 



1 Introduction 

Two-proton correlation functions at small relative mo- 
menta can probe the space-time extent of the reac- 
tion zone created in energetic heavy- ion reactions. This 
is due to the fact that the magnitude of nuclear and 
Coulomb final-state interactions (FSI) as well as anti- 
symmetrization effects depend on the spatial separation 
of the two protons during the emission process PEE 

Usu- 



ally, the correlation functions are evaluated as a func- 
tion of the magnitude q of the relative momentum vec- 



tor q — (pi — P2)/2. The interplay of the attractive S- 
wave nuclear interaction and the Coulomb repulsion and 
anti-symmetrization produce a minimum at q=0 and a 
maximum in the correlation function at g '--^ 20 MeV/c 
p. Most analyses give only upper limits for the spatial 
extent of the source due to the ambiguity of radius and 
lifetime of the source, i.e. model calculations simulating 
large sources with short lifetimes will give very similar cor- 
relation functions as model calculations simulating small 
sources with long lifetimes QIE]. 
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At beam energies below about 100 A-MeV, heavy-ion 
experiments performed mainly at the National Supercon- 
ducting Cyclotron Laboratory (NSCL) at Michigan State 

systematic investigations of pp small-angle correlations. 
Furthermore, new methods allowing to deduce the emis- 
sion source function from two-particle correlations have 
been proposed. Thus, recently, the technique of source 
imaging j23ll24ll25ll26| . i.e. the numerical inversion of the 
correlation function, has been applied successfully to pp 
correlations studied in heavy-ion experiments not only in 
the MSU-NSCL energy range 23 ,27 ,28^ but also at signif- 
icantly higher beam energies provided by the Brookhaven- 
AGS EH or even by the CERN-SPS 

However, in the energy domain of the heavy-ion syn- 
chrotron SIS at GSI Darmstadt, ranging from about 100 
to 2000 A-MeV, only few data on pp correlations are avail- 
able: proton-proton correlations at small relative momenta 
were measured with the "Plastic Ball" at BEVALAC for 
two systems, Ca on Ca and Nb on Nb, at 400 A-MeV 
as a function of proton multiplicity j^. Freeze-out densi- 
ties of about 25% of normal nuclear-matter density were 
extracted. Two-proton small-angle correlations were stud- 
ied at SATURNE in Ar -I- Au reactions at 200 A-MeV [TUI 
and in Ne- and Ar-nucleus collisions between 200 and 1000 
A-MeV Breakup densities and time scales in specta- 
tor fragmentation following the reaction Au -I- Au at 1 
A-GeV were investigated with the ALADiN spectrome- 
ter at SIS/GSI 18,19'. Finally, the space-time extent of 
the proton emitting source was studied with the An detec- 
tor FOPI measuring central collisions of Ni -I- Ni at 1930 
A-MeV and of Ru + Ru at 400 A-MeV IHIES- Conse- 
quently, the extension of the pp-correlation data base is 
highly desirable. 

Since the theoretical pp correlation function only re- 
quires the knowledge of the one-body phase-space distri- 
bution of protons, it is possible to generate small- angle 
correlations with any microscopic theory that calculates 
the time evolution of the one-body distribution function 
In addition to the usage of Gaussian source distri- 
butions we will also follow this procedure since it allows 
a more adequate comparison of experimental and theo- 
retical correlation functions. Microscopic transport mod- 
els, hke Boltzmann-Uehling-Uhlenbeck (BUU) piTllg] or 
Quantum Molecular Dynamics (QMD) [31] in their vari- 
ous implementations, do not only allow to trace the ex- 
pectation values of the positions of the various particles 
in coordinate space but also the corresponding time in- 
formation. Moreover, they carry the correlation of coordi- 
nate and momentum space coordinates of the individual 
particles. This correlation is generated intrinsically dur- 
ing the expansion following the compression phase of the 
heavy-ion collision. On a macroscopic scale, the expan- 
sion constitutes a collective motion of the various particle 
species, called radial flow. Thus, it is natural to expect 
that a hybrid model consisting of BUU approach plus FSI 
model should be able to reproduce the experimental data 
and give the relevant insight into the physics of the heavy- 
ion collision until particle freeze out. Furthermore, there is 



reasonable hope that a systematic comparison of as much 
as possible experimental information, e.g. on sideward j32| 
and elliptic flow on nuclear stopping 34_ and cluster 
production |35j . on the shape (in momentum space) of 
the participant source |36| or on charged pion production 
|37| . with transport model predictions will allow to con- 
strain the essential input parameters of the model, e.g. the 
equation of state (EoS) and/or the nucleon-nucleon cross 
section within the nuclear medium. Earlier it was reported 
that the pp correlation function generated for ^''N-l-^^Al 
collisions at 75 A-MeV is only weakly dependent on the 
EoS OEj- It is not clear whether the sensitivity to the 
microscopic details of the transport approach increases 
for higher beam energies and/or heavier projectile-target 
combinations. Therefore, in the following we will study 
in a systematic way how the BUU+FSI hybrid model re- 
produces the data and what kind of information can be 
gained. 

The paper is structured as follows: In sect.El the ex- 
perimental basis is described shortly, the involved event 
classes are explained and the correlation function is de- 
fined. In sect.|3|the experimental pp correlation functions 
with their dependences on system size and beam energy 
are presented and compared to model predictions. Finally, 
the results are summarized in sect.0] 



2 The experiment 

The experiments have been performed at the heavy-ion 
synchrotron SIS at GSI Darmstadt. Symmetric collisions 
are carried out by irradiating targets of 1 % interaction 
thickness of "^^Ca and ^^^Au with the corresponding ions 
of beam energies between 400 and 1500 A-MeV. At a pro- 
jectile energy of 400 A-MeV the system of intermediate 
mass, ^^Ru -I- ^^Ru, is included which was explored ex- 
tensively with respect to the space-time difference of light 
charged particles emission in central collisions |22| . 

The present analysis uses a subsample of the data, 
taken with the outer Plastic Wall/Hclitron combination 
of the FOPI detector system The Plastic Wall 

delivers - via energy loss vs. time-of-flight (TOF) mea- 
surement - the nuclear charge Z and the velocity /3 of the 
produced particles and the corresponding hit positions. 
The Helitron gives the curvature (which is a measure of 
momentum over charge (p/Z)) of the particle track in the 
field of a large super-conducting solenoid. Since the mo- 
mentum resolution of the Helitron is rather moderate, this 
detector component serves for particle identification only. 
The mass m is determined via mc — (p/ Z)Hei /{Pl/Z)pia, 
where 7 = (l-/?^)-^/^. The Plastic WaU and the Helitron 
have full overlap only for polar angles between 8.5 degrees 
and 26.5 degrees. The corresponding flight paths amount 
to 450 cm and 380 cm, respectively. Monte-Carlo simula- 
tions have been performed in order to study the influence 
of the finite detector granularity and of the TOF and po- 
sition resolutions on the velocity and finally on the proton 
momentum. The resolution of both quantities is governed 
by the TOF resolution, which is (Jtof = 80 (120) ps for 
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short (long) scintillator strips located at small (large) po- 
lar angles |3H1- The velocity can be determined with a 
precision oi 6$/^ < 1 %. 



2.1 Event classification 

Typically, a few hundred thousand to a million central 
events are collected for each individual collision system 
by demanding large charged-particle multiplicities to be 
measured in the outer Plastic Wall. The corresponding 
integrated cross sections for the collision systems Ca + Ca, 
Ru + Ru, and Au + Au comprise about 15 %, 10 % and 
10 % of the total cross section, respectively. For the present 
central-event class one would expect - within a geometrical 
picture - average impact parameters of about 2 fm, 2 fm, 
and 3 fm. Taking into account typical dispersions of the 
impact parameter distribution of about 1 fm as found from 
GEANT simulations .4Q.!, the number of participants in 
the central source is estimated to Apart = 56 ± 12 (147 ± 
21, 302 ± 34) for the system Ca + Ca (Ru + Ru, Au 
+ Au). More peripheral collisions have been measured, 
too. But, due to both the strong downscaling (typically a 
factor of 16) of the corresponding online trigger and the 
lower proton multiplicities, the pp pair statistics was not 
sufficient for correlation analyses. 

In previous investigations of central Au + Au collisions 
between 100 and 400 ^-MeV beam energy it was found 
that the correlation function of pairs of intermediate mass 
fragments is strongly affected by the collective directed 
sideward flow of nuclear matter |41ll42j . This directed side- 
flow causes an enhancement of correlations at small rel- 
ative momenta. The enhancement results from mixing of 
differently azimuthally oriented events; it vanishes if the 
events are rotated into a unique reaction plane, which is 
determined by the standard transverse momentum anal- 
ysis 031 ■ Consequently, the technique of event rotation is 
applied also to the present data in order to prevent that 
such artificial correlations are introduced into the refer- 
ence momentum distribution of the correlation function 
(cf. sect.E21). 



2.2 Correlation function 

Let Yi2(pi,P2) be the coincidence yield of pairs of par- 
ticles having momenta pi and p2. Then the two-particle 
correlation function is defined as 



l-t-R(pi,P2) =AA 



events, pairs 



yi2(Pl,P2) 



E 



events,pairs 



.mix (Pl,P2) 



The sum runs over all events fulfilling the above mentioned 
global selection criterion and over all pairs satisfying cer- 
tain conditions given below. Event mixing, denoted by the 
subscript "mix", means to take particle 1 and particle 2 
from different events. The normalization factor J\f is fixed 
by the requirement to have the same number of true and 
mixed pairs. The statistical errors of all the correlation 
functions presented below are governed by those of the 



coincidence yield, since the mixed yield is generated with 
two orders of magnitude higher statistics. The correlation 
function is then projected onto the relative momentum 
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\pr~pT\ 



(2) 



Here, p^"* are the particle momenta calculated in the cm. 
system of the colliding nuclei. From the velocity resolution 
as estimated in sect.|21 the corresponding q resolution is 
deduced to be 5g = (4 ± 1) MeV/c. A similar value can 
be derived directly from the experimental dispersion of the 
resonance due to the narrow 2.186 MeV state of ^Li ( J'^ — 
3+, r' = 24 keV) in the a-d correlation function 22 . All 
theoretical correlation functions presented in sect.|31 are 
folded with this q resolution. 

As in previous proton-proton correlation analyses |2()l 
I21ll22j an enhanced coincidence yield at very small rela- 
tive angles is observed, which is due to double counting 
caused mainly by scattering in the scintillator strips. This 
disturbing yield is reduced strongly by the requirement to 
match the particle hits on the Plastic Wall with the cor- 
responding tracks in the forward drift chamber Helitron. 
The remaining left-over of doubly counted scattered par- 
ticles is eliminated by excluding, around a given hit, posi- 
tions within a rectangular segment of azimuthal and po- 
lar angle differences - 02| < 4° and \di ~ < 2°. 
The same procedure is applied to the uncorrelated back- 
ground, hence keeping the influence of the exclusion onto 
the correlation function as small as possible. GEANT sim- 
ulations 40 have shown that at very small relative mo- 
menta, q < (12 — 15) MeV/c, still a small bias of the 
correlation function can not be excluded |20II22| . Thus, 
the corresponding regions in the correlation functions are 
not taken into consideration when comparing the experi- 
mental data with model predictions. 



3 Results and comparison with model 
predictions 

For description of the FSI model we refer to ref. and 
references cited therein. As source distributions we use ei- 
ther Gaussian density profiles in coordinate and momen- 
tum space or the output of the BUU transport model. The 
theoretical correlation function is given as 



l + i?(P,q) 



d^rS{r,P)\^^{r) 



(3) 



Here, the wave function tZ^ describes the relative motion 
of the two emitted particles and S is the source function. 
Gaussian sources are defined in [22] • In case of calculat- 
ing the single-particle phase-space distribution within the 
BUU model the source function is given as 

5(r,P)=AAp^<5(|-p05(^-p,)x 

hi 

S{r-{r,-r,) + ^{U~t,)). (4) 
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Ti, ti, and Pi are the space, time and momentum coor- 
dinates of particle i at freeze out. The parameter P rep- 
resents the sum momentum of the observed particle pair. 
The normalization J\fp is chosen such that 

J d^rS{r,P) = 1. (5) 

To allow the summation over a sufficient number of par- 
ticles, the 6 function treating the momentum is replaced 
by a Gaussian with dispersion Ap, 

Note that all coordinates are calculated in the cm. sys- 
tem of the colliding nuclei. The relevant BUU input pa- 
rameters are i) the number of test particles per nucleon 
{i.e. parallel runs or pseudo events, here chosen to 200), 
ii) nuclear charge and mass numbers of target and pro- 
jectile, At, Zt, Ap, Zp, iii) the projectile energy Eproj, 
iv) the impact parameter b, and v) the stiffness of the 
Equation of State (EoS), either hard (k = 380 MeV), 
medium {k = 290 MeV) or soft {k = 215 MeV) 44 . 
A proton is considered "emitted" when the surrounding 
density falls below a value of Pcutoff = 0.02/fm'^ and 
when subsequent interaction with the mean field does not 
cause recapture into regions of higher density until the 
calculation is terminated at 150 fm/c. The cutoff value of 
about one eighth of normal nuclear matter density implies 
ceasing NN interactions in a sufficiently diluted nuclear 
matter. This choice was approved in systematic inves- 
tigations of proton-proton correlations in ^''N-l-^^Al and 
i^N+i^^Au reactions at 75 A-MeY 0, in ^GAr+^^Sc col- 
hsions at 80 A-MeV fTTltn] and in "^^Ar+i^^Au reactions 
at 200 ^-MeV [E]. 

Finally, for each run and each proton, the sep- 
tuplet of space-time-momentum coordinates Ni = 
{x,y, z,t,Px,Py,Pz)i is taken from the BUU approach at 
the corresponding freeze-out time and then processed 
within the FSI model. The relevant parameters of the 
FSI model are the so-called "observation quantities", i.e. 
the sum momentum P, the momentum dispersion Ap 
and the polar angle 9™^. These quantities are taken di- 
rectly from the experimental distributions in the cm. sys- 
tem, e.g. P = 2(p=™). Typical values are O""" = 45°, 
Ap = 150 MeV/c, and (p"") = 300 (400) MeV/c for 400 
(1500) j4-MeV beam energy. 

3.1 System-size dependence 

Figure^ shows, for 400 ^-MeV beam energy, the system- 
size dependence of the pp correlation function. Clearly, the 
pp-correlation peak decreases with increasing system size 
and almost vanishes for the largest collision system Au 
-|- Au. The peaks {i.e. data within a relative-momentum 
range of 15 < g/MeV/c< 50) of the experimental cor- 
relation functions (symbols) are best described by the 
theoretical ones (lines) for apparent Gaussian radii of 
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Fig. 1. Correlation functions of proton pairs from central 
collisions of Ca -I- Ca, Ru -I- Ru, and Au + Au at 400 A-MeV. 
Experimental data (symbols) are compared to predictions of 
the FSI model with Gaussian sources and zero lifetime (lines). 
The corresponding apparent source radii are indicated. 

^ (3.0 ± 0.1) fm for Ca + Ca, (4.0 ± 0.15) fm for Ru 
+ Ru and (5.1 ± 0.2) fm for Au + Au colhsions. Note 
that the term "apparent" means that neither the con- 
tribution of the finite emission duration (inappropriately 
called "lifetime", leading to an increase of the effective 
source radius) nor the effect of collective radial expansion 
of the participant zone (leading - via strong correlations 
of coordinate and momentum space - to a reduction of 
the effective source size) are taken into account [201 • In 
the following, the asterisk indicates apparent quantities, 
e.g. the sharp-sphere and r.m.s. radii i?*^ = Rq and 
^rms = VSRq, respectively. For Gaussian source distri- 
butions the effective radius can be written as 



^o = \/tt; + («^)'' (7) 

where t is the duration of emission, v — P/2m is the pair 
velocity, and e — ejiow/^therm represents the ratio of radial 
flow energy e/ioiu and the energy of the random thermal 
motion etherm = ^T. 

For both the small and the medium-size systems, our 
results can be compared directly with existing data: Gaus- 
sian radii rg were derived from pp correlations measured 
with the "Plastic Ball" at BEVALAC for the reactions 
Ca -t- Ca and Nb + Nb, both at 400 MeV per nucleon, as 
function of proton multiplicity (H]. Taking into account the 
different radius definitions, i.e. Vg = \^Rq, the apparent 
source radii for central collisions agree perfectly. Conse- 
quently, the freeze-out density of about 25% of normal 
nuclear matter density deduced in ref. 9 is in accordance 
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400 AMeV A+A, BUU coordinate space 
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Fig. 2. Coordinate-space distributions of freeze-out points 
of protons from central collisions of Ca -I- Ca (top), Ru + 
Ru (center), and Au 4- Au (bottom) at 400 ^-MeV simulated 
by the BUU transport approach. The corresponding impact 
parameters are indicated at the right margin. The left, middle 
and right columns give the distributions in the (y,x), {x,z) 
and (y, z) plane, respectively. 

with our results presented in sect . 13. 21 below. Instead of 
disentangling the space-time-flow ambiguity of eq. ^20 
I22j , in the following we try to reproduce the experimental 
correlation functions by BUU-I-FSI model calculations. 

Before we confront the data with theoretical correla- 
tion functions we briefly focus onto the coordinate-space 
distributions of freeze-out points from BUU calculations 
displayed in fig.El for the three systems at 400 A-MeV 
beam energy. Note that the coordinates are not taken at 
fixed time but when the local density drops below the cut- 
off value. As expected, the proton source increases with in- 
creasing system size. Additionally to the coordinate-space 
distribution shown in fig.|3 as an example for the Au -f 
Au system fig-El displays the momentum distribution pro- 
jected onto the reaction plane (j>x,Pz), the {y,x) distri- 
bution for a central z slice and the radial density profile. 
Here, the effective density is given in units of normal nu- 
clear matter density, po = 0.168/fm^. The corresponding 
r.m.s. radius Rrms = a/ (r^) of the proton emitting source 
is indicated by an arrow. No essential changes are observed 
from these phase-space plots when increasing the beam en- 
ergy, except a trivial extension in momentum space. Since 
even the r.m.s. radius of the breakup configuration (cf. ar- 
row in lower right panel) is almost the same for 400 and 
1500 j4-MeV beam energy, one would expect a very sim- 
ilar height of the pp correlation peaks. However, also the 
correlation of coordinate and momentum space, constitut- 
ing a fast collective motion known as radial expansion, is 
expected to increase with beam energy. It should lead to 



Fig. 3. Phase-space distribution of freeze-out points of pro- 
tons from central collisions of Au -I- Au at 400 A-MeV simu- 
lated by the BUU transport approach for an impact parameter 
of & = 3 fm. Left column (from top to bottom): The distribu- 
tions in the {px,Pz), {x,z), and {y,z) plane. Right column: 
The {y,x) distribution in total (top) and for a central z slice, 
|2:| < 2 fm (center). The bottom panel gives the corresponding 
radial density profile (in units of the saturation density). The 
arrow indicates the r.m.s. radius of the source. 

a reduction of the apparent source radius and hence to an 
increase of the pp peak height (cf. cq. O). This presump- 
tion will be verified in sect. l3.2l 

For 400 ^-MeV beam energy, fig.^ compares the ex- 
perimental system-size dependence of the pp correlation 
function with the corresponding output of the BUU-I-FSI 
model. Hardly any differences have been found between 
the pp correlation functions derived from BUU simula- 
tions using the three different parameterizations of the 
EoS as can be inferred from the narrow hatched bands in 
Fig-El (Similar findings have been reported for ^"'N-l-^^Al 
collisions at 75 A-MeV |3|H1 and for '^^Ar-H^^^Au reactions 
at 200 A-MeV UU].) Therefore, in the following we wiU re- 
strict ourselves to the medium stiff EoS. 



3.2 Beam-energy dependence 

Figure|Sl shows the excitation function of experimental pp 
correlations from central collisions of Ca + Ca at beam 
energies from 400 to 1500 A-MeV (symbols). A strong in- 
crease of the correlation peak is obvious corresponding 
to a decrease of the apparent source size. Indeed, when 
comparing with the theoretical correlation functions de- 
rived from the FSI model with static Gaussian sources 
and zero lifetime (lines), the experimental pp peak is best 
reproduced with apparent radii systematically decreasing 
from 3.0 fm at 400 A-MeV to 2.5 fm at 1500 A-MeV. This 
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Fig. 4. The same as fig.0 but here the experimental corre- 
lation functions (symbols) are compared to predictions of the 
BUU+FSI model (lines). The narrow hatched bands around 
the lines indicate the systematic variation of the calculations 
for the three different parameterizations of the EoS (see text). 



Ca+Ca, central collisions 



2 




0.02 0.04 0.06 0.08 0.1 0.12 



q/(GeV/c) 

Fig. 5. Correlation functions of proton pairs from central 
collisions of Ca + Ca at projectile energies from 400 to 1500 
A-MeV. Experimental data (symbols) are compared to predic- 
tions of the FSI model with Gaussian source and zero lifetime 
(lines) . The corresponding apparent source radii are indicated. 



shrinking of the apparent source has to be traced back 
to two facts. First of all, at higher beam energies when 
the heavy-ion coUisions is processing faster, the contri- 
bution of the emission duration onto the effective source 
radius is smaller. Usually, finite emission times imply a 
reduced Pauli suppression in emission direction and hence 
lead to an increase of the effective source radius 1.201.) 
Typically, for the energy range we are dealing with in 
the present investigation, these times are in the order of 
10 fm/c pnil2m22| . Secondly, also the collective (radial) 
expansion increases with increasing projectile energies. It 
is well known that strong correlations of coordinate and 
momentum space lead to a drastic reduction of the appar- 
ent source radius as derived from a comparison of the data 
with results of a FSI model incorporating the Koonin- 
Pratt formalism with static Gaussian sources and finite 
life (emission) times. 30 % smaller radii are observed |20l 
I21II22| for a ratio of collective to random thermal energies 
of unity ^5] . 

FigureEl shows the excitation function of experimen- 
tal pp correlations from central collisions of Au -I- Au 
(symbols). From comparisons with the output of the FSI 
model (lines), we find an optimum reproduction of the 
experimental data in the peak region for apparent radii 
dropping from 5.1 fm at 400 A-MeV to 4.0 fm at 1500 
A-MeV. The trend found in central Au -I- Au collisions 
is the same as for the smaller system Ca -I- Ca, but with 
lower peak values. This effect is expected due to the much 
larger number of involved nucleons in central Au -I- Au 
collisions. Looking at the volumes V* = 4tt{RI^)^ /3 de- 
rived from the Gaussian radii given in figs.m3and|ni our 



expectations are confirmed: for the same beam energy, the 
volume ratio is approximately equal to the system size ra- 
tio of 197/40. Our investigation now points to the breakup 
densities derived from the given radii. With the given vol- 
umes and the number of participants estimated in sect. 12. ll 
apparent breakup densities (normalized to the density 
of normal nuclear matter) of p* / po = (j^part/V*)/ po — 
0.27±0.06 (0.29±0.05, 0.29±0.05) are deduced for Ca -)- 
Ca (Ru + Ru, Au + Au) collisions at 400 ^-MeV. These 
numbers increase to values of 0.47 ± 0.11 (0.60 ± 0.10) 
for Ca -I- Ca (Au + Au) at beam energies of 1500 A-MeV. 
Here, we remind once more to the fact that these apparent 
breakup densities are not corrected for both the effects of 
collective radial expansion and finite duration of emission 
(cf. eq. 0). 

Finally, we confront the experimental beam-energy de- 
pendence of two-proton correlations with the predictions 
of the BUU-I-FSI hybrid model. Figure|71 shows the ex- 
citation function of pp correlations for central Ca + Ca 
collisions in comparison to the model results. With the 
standard parameters in the transport approach the ex- 
citation functions are surprisingly well reproduced. Also 
for the larger system, Au -I- Au, the energy dependence is 
nicely described (cf. fig.|Hl). Here, the experimental pp cor- 
relation peaks are systematically slightly underestimated. 

We want to point out that the agreement between our 
experimental correlation functions and those predicted by 
BUU is much better than in previous experiments per- 
formed at lower beam energies. Thus, pp correlation func- 
tions predicted by BUU calculations overestimated the 
measured central collision data of the reactions Ar + Au 



7 



Au+Au, central collisions 




0.02 0.04 0.06 0.08 0.1 0.12 



q/(GeV/c) 

Fig. 6. The same as figure|3 but for Au + Au collisions. 

at 200 A-MeV (TUI and of Ar + Sc at 120 and 160 A-MeV 
|16| . This deficiency of the BUU model was attributed to 
its inability to treat the population of particle unbound 
resonances and their decay via delayed particle emission 
|16II28| . Obviously, at our higher beam energies this slow- 
emission component of the source function becomes unim- 
portant, presumably since the number of heavy fragments 
drops strongly with energy |45| . Thus, the fast expansion- 
explosion scenario seems to be well modelled within the 
transport approach. 

4 Summary 

In the present paper, small angle correlations of pairs of 
protons emitted in central heavy-ion collisions at beam 
energies from 400 to 1500 ^-MeV are investigated. New 
data are presented which substantially enlarge the exper- 
imental data basis on pp correlations in the SIS energy 
range. 

The system-size and beam-energy dependences of the 
pp correlation peak are studied. The peak decreases, i.e. 
the apparent source radius increases, with increasing sys- 
tem size. This system-size dependence is expected from 
the larger volume due to a larger number of participants 
in the central source. With increasing beam energy the 
peak of the pp correlation function increases, and hence 
the apparent source radius decreases. This behaviour is 
interpreted as the common action of shorter time scales 
and stronger collective radial expansion. Correlations of 
coordinate and momentum space constituting the latter 
radial flow are known to be responsible for the apparent 
shrinking of the source size found in small-angle correla- 
tion studies. 

Both findings, shorter emission durations due to 
shorter time scales of the central heavy-ion collision and 



Ca+Ca, central collisions 
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Fig. 7. Correlation functions of proton pairs from central 
collisions of Ca + Ca at projectile energies from 400 to 1500 
A-MeV. Experimental data (symbols) are compared to predic- 
tions of the BUU+FSI model (lines). 

stronger radial expansion due to - increasing with pres- 
sure - space-momentum correlations, are well incorpo- 
rated within microscopic transport models. Accordingly, 
using the Boltzmann-Uehling-Uhlenbeck transport ap- 
proach, the phase-space points at freeze out are calculated 
and afterwards processed within a final-state-interaction 
model providing the two-proton correlation function. This 
BUU+FSI model well reproduces the system-size and 
beam-energy dependences of the experimental correlation 
function. However, the pp correlation function is found 
rather insensitive to the stiffness of the equation of state 
used in the transport model. 

We are grateful for many discussions to H.W. Barz, B. Kampfer 
and F. Dohrmann. 
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